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MODELING  BACKSCATTERING  FROM  A  ROUGH  SEAFLOOR 
WITH  SEDIMENT  INHOMOGENEITIES 


Jorge  C.  Novarini,  Planning  Systems,  Inc.,  Slidell,  LA  70458 
Jerald  W.  Caruthers,  Naval  Research  Laboratory, 
Stennis  Space  Center,  MS  39529 


ABSTRACT 

Current  models  used  to  predict  the  backscattering  strength  of  the  ocean  floor  are 
either  very  involved,  requiring  geoacoustic  parameters  usually  unavailable  for  the  site 
in  practical  applications,  or  overly  simplistic,  relying  mainly  on  empirical  terms  such  as 
Lambert’s  law.  In  any  case,  solutions  are  very  approximate  and  the  problem  is  still  far 
from  being  solved.  In  this  work,  a  model  is  presented  that  avoids  empirical  functional 
forms,  yet  requires  only  a  few  physical  parameters  to  describe  the  surficial  sediments, 
often  tabulated  for  typical  sediments.  The  aim  of  the  work  is  to  develop  a  simple 
algorithm  for  operational  prediction  of  bottom  reverberation  with  only  one  free 
parameter,  i.e.,  the  volume  scattering  coefficient.  The  algorithm  combines  a  two-scale 
surface  scattering  model  [Caruthers  and  Novarini,  IEEE  J.  Oceanic  Engr,  v  18,  p  100- 
107,  1993  and  Novarini  and  Caruthers,  IEEE,  J.  Oceanic  Engr,  v  19,  p  200-207,  1994], 
with  scattered  contributions  originating  from  inhomogeneities  within  the  sediments, 
taking  into  consideration  the  rough  interface  [Ivakin  and  Lysanov,  Sov.Phys.  Acoust. 
v  27,  p  212-215,  1981].  No  specific  mechanism  is  assumed  for  scattering  at  the 
volume  inhomogeneities;  however,  the  inhomogeneities  are  assumed  to  be  uniform 
and  isotropic.  The  volume  scattering  coefficient,  combined  with  the  bottom  attenuation 
and  density  and  referenced  to  the  surface,  plays  a  role  similar  to  the  Lamberts 
constant  in  empirical  models.  The  model  is  exercised  on  a  variety  of  published 
datasets  for  low  and  moderately  high  frequency.  In  general,  the  model  performs  very 
well  for  both  fast  and  slow  sediments,  showing  a  definite  improvement  over  Lambert  s 

law. 


I.  INTRODUCTION 


Scattering  from  the  seafloor  has  been  the  subject  of  many  papers  in  the  last 
forty  years.  However,  in  spite  of  these  many  efforts,  the  problem  is  far  from  being 
solved  Acoustic  bottom  interaction  is  an  extremely  complicated  physical  phenomenon, 
since  it  involves  scattering  at  a  statistically  rough  penetrable  surface  overlaying  an 
inhomogeneous  viscoelastic  medium.  Each  of  these  phenomena  are  individually 
complex  and  their  combined  effect  is  almost  unmanageable.  An  additional 
complication  is  that  some  characteristics  of  the  sediments  are  difficult  to  measure,  and 
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detailed  bottom  structure  is  known  only  over  very  limited  geographical  regions. 
Moreover  acoustic  experiments  on  bottom  interaction  are  not,  as  a  rule,  accompame 
by  the  geophysical  information  necessary  to  model  adequately  the  data  or  interpret 
the  results.  As  a  consequence,  models  may  quickly  become  complicated  and 

justification  weak. 

In  the  past  bottom  reverberation  was  considered  to  be  produced  mainly  at  the 
water-sediment  interface,  with  possible  contributions  from  subbottom  rough  interfaces 
In  the  last  few  years  evidence  has  accumulated  which  indicates  that  scattering  a 
volume  inhomogeneities  within  the  sediment  is  a  strong  contributor,  and  this  has 
stimulated  an  increased  interest  in  the  subject.1'6  Although  early  attempts  were  made 
to  evaluate  the  scattering  strength  produced  within  the  sediment,  these  dealt  with 
volume  scattering  from  a  medium  bounded  by  a  flat  interface  (Stockhausen  and 
Zhitkovskii8).  In  a  clever  and  elegant  paper,  Crowther  considered  both  bottom 
roughness  and  volume  scattering;  but  when  dealing  with  the  latter  contribution  a  flat 

interface  was  implicitly  assumed. 

The  pioneer  work  on  the  volume  scattering  contribution,  including  the  effects 
of  the  rough  boundary,  can  be  credited  to  Ivakin  and  Lysanov.  Their  volume 
scattering  component  was  based  on  an  isotropic  distribution  of  scatters  beneath  a 
surface  that  was  rough  in  the  large  scale.  The  contribution  of  volume  scatters  was 
averaged  over  the  distribution  of  slopes.  The  averaging  process  was  a  little  awkward 
and  they  did  not  include  a  separate  surface  roughness  scattering  mechanism. 
Following  the  approach  of  Ivakin  and  Lysanov,  Jackson  and  his  co-workers 
developed  a  model  that  included  separate  surface  scattering  terms.  They  adopted 
Stockhausen’s  model  for  the  volume  contribution,  and  carried  out  the  large-scale  slope 
averaging  process  numerically.  Recently,  Lyon  et  al.  extended  Jackson  s  work  to 
include  specific  volume  scattering  mechanisms  and  scattering  at  the  subbottom 

interface. 

Most  current  efforts  in  the  research  community  are  concentrated  on  modeling 
bottom  backscattering  making  use  of  an  increasing  number  of  geoacoustic  parameters 
(e.g.,  sound  speed  gradient  in  the  sediments,  fluctuations  in  compressibility  and  in  the 
density,  stratification,  attenuation).  However,  most  of  those  parameters  are  usually 
unavailable  and,  as  a  matter  of  fact,  no  firm  experimental  evidence  exists  to  support 
any  specific  scattering  mechanism  requiring  such  a  detailed  description. 

The  intricacies  of  this  approach  stand  in  stark  contrast  to  the  simplicity  of 
Lambert’s  law,  which  is  currently  used  for  most  sonar  applications.  Lamberts  law 
requires  a  single  empirical  coefficient,  which  was  established  by  Mackenzie,  based 
on  a  limited  dataset  taken  in  the  early  1960’s,  to  be  -27  dB.  Work  is  currently  in 
progress  to  develop  values  for  the  Lambert  constant  on  a  regional  basis  (Geddes  ). 
However,  Lambert’s  law  has  shortcomings.  It  has  a  cosine-squared  angular  functional 
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form  which  sometimes  fails  to  fit  data,  especially  at  lower  grazing  angles,  and  it  is 
inappropriate  for  near-specular  scatter  (i.e.,  near-normal  backscatter  and  forward 
scatter).  To  mediate  against  this  latter  problem,  Lambert’s  law  is  sometimes  combined 
with  some  near-specular  scattering  algorithm.13,14  One  of  these  algorithms,  the  rough 
facet  model  (RFM),  uses  a  two-scale  surface  scattering  model  for  predicting  the 
specular  lobe,  and  resorts  to  Lambert’s  law  for  directions  away  from  specular. 

To  replace  the  Lambert  term,  this  present  work  combines  the  rough  facet  model 
with  a  volume  scattering  component  based  on  the  work  of  Ivakin  and  Lysanov,  and 
implements  resonant  (Bragg)  scattering  at  the  microroughness  scale.  The  seafloor  is 
assumed  to  be  a  fluid  containing  isotropic  volume  inhomogeneities  witn  no 
stratification.  The  model  uses  the  volume  scattering  cross-section  of  the  sediment  as 
a  free  parameter  which,  combined  with  the  sediment  attenuation  and  density,  plays  a 
role  analogous  to  the  Lambert  constant.  The  basic  aim  of  this  work  is  to  provide  a 
simple  model  for  predicting  the  bottom  backscattering  strength  which  avoids  empirical 
expressions  (such  as  Lambert’s  law)  while  keeping  the  calculations  and  input 
requirements  at  a  minimum.  No  attempt  is  made  here  to  connect  the  resulting  volume 
scattering  strength  with  a  specific  volume  scattering  mechanism.  It  is  shown  that  the 
overall  performance  is  satisfactory,  and  constitutes  a  definite  improvement  over 
Lambert’s  law  in  fitting  a  variety  of  datasets. 


II.  DESCRIPTION  OF  THE  SCATTERING  MODEL 

The  backscattering  coefficient  of  the  seafloor,  mbs,  is  assumed  to  be  made-up 
of  three  contributions: 

mbs  =  mf  +  m^  +  mv  ,  (1) 

where  mf  represents  the  scattering  contribution  from  fine-scale  components  of  the 
bottom  roughness  (geomorphology  at  a  scale  just  below  large-scale,  deterministic 
bathymetry  but  larger  than  an  acoustic  wavelength)  within  the  footprint  of 
ensonification;  m  is  the  contribution  from  the  micro-scale  roughness 
(microroughness),  and  mv  is  the  contribution  from  the  inhomogeneities  within  the 
volume  of  the  sediments,  with  due  regard  to  the  fact  that  it  comes  through  the  rough 
interface.  Henceforth,  this  last  contribution  will  be  referred  to  as  the  "volume-induced 
bottom  scattering  coefficient".  The  large-scale  roughness  components,  associated  with 
long  surface  wavelengths  present  in  bathymetric  databases  are  handled 
deterministically  and  are  excluded  from  this  analysis  (i.e.,  as  atilt  in  the  footprint  which 
is  removed  deterministically).  The  quantity  10log(mbs)  is  the  bottom  backscattering 
strength  in  dB.  The  individual  terms  in  Eq.  (1)  may  also  be  expressed  in  dB.  Generally, 
there  is  only  a  small  angular  range  in  which  the  terms  must  be  added  in  their 
magnitude  quantities. 
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Thfi  Rough  Facet  Mode] 


The  first  term  (mf)  is  handled  by  the  rough  facet  model.  This  term  represents 
the  collective  contribution  from  reflections  at  the  fine-scale  facets  containing 
microroughness.  That  is,  it  is  controlled  by  a  facet  reflection  process  which  is  partially 
coherent  for  individual  facets  but  collectively  incoherent  when  summed  over  the  sonar 
footprint  Additional  fundamental  assumptions  in  approach  are  that  the  facet  slopes 
have  a  gaussian  distribution  over  the  footprint  and  that  the  source  and  receiver  are  in 
the  farfield  of  the  footprint.15  The  term  has  the  form 


mf  =  M(k,a  6)  •  F(6fl0)  , 


(2) 


where  M(-)  is  the  coherent  reflection  coefficient  due  to  the  microroughness  on  a  facet 
and  F(-)  is  the  contribution  of  the  fine-scale  slopes  of  the  facets.  The  additional 
parameters  are:  k  is  the  acoustic  wavenumber;  is  the  rms  height  of  the 
microroughness;  6f  is  the  rms  slope  of  the  fine-scale  facets;  and  e  is  the  local  angle 
of  incidence  measured  from  normal.  (We  will  later  have  use  for  the  rms  slope  angle 
defined  by  af=tan'16f.)  The  coherent  reflection  coefficient  is  given  by 


M  =  exp(-g  J  ,  W 

where  g^  is  the  roughness  parameter  of  the  microroughness,  i.e., 

=  4  ct^2  k2  cos20  .  (4) 

The  function  F(-)  is  the  high-frequency  limit  of  the  Helmholtz/Kirchhoff  theory.  5  For 
the  case  of  monostatic  backscattering  this  reduces  to 


F  = 


(5) 


where  R  is  the  pressure-reflection  coefficient  of  the  sediment,  assuming  it  to  be  a 
fluid.  As 'in  any  two-scale  roughness  theory,  the  choice  of  the  partition  wavenumber 
between  the  fine-scale  roughness  and  the  microroughness  is  crucial  to  the  results. 

Power  spectra  used  to  describe  seafloor  roughness  are  most  often  of  the 
isotropic  power-law  form, 


where  K  is  the  spatial  frequency  of  the  surface  roughness  (i.e.,  the  reciprocal 
wavelength  in  cycles/meter).  (Note  our  use  of  the  spatial  frequency  K  for  surface 
specifications  and  the  wavenumber  k  for  acoustics  applications.) 

It  should  be  noted  that,  in  the  rough  facet  model,  the  two  surface  parameters 
(the  rms  height  of  the  micro-scale  surface  and  the  rms  slope  of  the  fine-scale  surface) 
are  bandlimited  quantities.  The  variance  of  the  microroughness  for  a  2-D  isotropic 

surface  is  given  by 


2  _ 


=  2n 


W(K)  K  dK 


(7) 


The  variance  of  the  slopes  of  the  fine-scale  surface  by 


6f  =  2rr 


K, 


K, 


W(K)  K3  dK  . 


(8) 


In  Eqs.  (7)  and  (8),  K,  and  Kh  indicate  the  low  and  high  spatial  frequencies  present  on 
the  surface,  respectively.  These  spatial  frequencies  are  determined  by  the  size  of  the 
footprint  and  the  grid  spacing,  respectively.  Kc  denotes  the  partition  (cut-off)  spatial 
frequency.  For  analytical  work  Kh  can  approach  infinity;  but  a  K,  approaching  zero 
would  give  rise  to  problems  with  the  power-law  spectrum,  and  is  not  physical  because 
the  footprint  of  a  sonar  is  always  limited  in  size. 

A  recent  numerical  study16  has  shown  that  the  appropriate  partition  for  the 
rough  facet  model  is  that  which  makes  g^  equal  to  unity.  Under  that  condition,  RFM 
models  the  specular  lobe  with  high  accuracy. 

The  reflection  coefficient  is  given  by 

R0  =  (u-1)/(u+1)  ,  (9) 


where 


u  =  m  cose  /  (n2  -  sin2e)Vz . 


(10) 


Since  the  sediment  is  considered  a  lossy  medium  in  this  work,  the  index  of  refraction 
n  is  complex,  i.e., 
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n  =  c/c1  =  n0(1+i^>),  and  m  =  p1  /  p  , 

where  c  is  velocity  and  p  is  density  (parameters  subscripted  with  "1"  here  and 
throughout  are  for  the  sediments).  The  parameter  <p  is  the  so-called  "loss  tangent" 
given  by 


=  Kp  log10(c1)  /  40 tt  , 

where  k  is  the  attenuation  factor  for  compressional  waves,  defined  as  the  attenuation 
coefficient  in  the  sediment  (dB/m)  divided  by  frequency.17  It  has  the  units  dB/(m-kHz) 
and  frequency  expressed  in  kHz. 

Microrouahness  Scattering 

In  Eq.  (2),  the  microroughness  is  only  present  as  a  loss  mechanism  in  the 
specular  direction  (through  the  factor  M).  However,  the  microroughness  re-scatters 
energy  in  all  directions,  including  the  receiver  direction,  and  its  contribution  to  the  total 
scattering  coefficient  must  be  accounted  for.  Since,  in  the  micro-scale,  we  are  dealing 
with  an  acoustically  smooth  surface  (ka^  <  1),  the  contribution  can  be  evaluated 
through  perturbation  theory,  leading  to  resonant  (Bragg)  scattering.  For  this  work  we 
have  adopted  the  expression  given  by  Brekhovskikh  and  Lysanov  for  Bragg 
scattering,  modified  for  our  definition  of  variance  (Eq.  (7)).  It  is  less  desirable  to 
assume  the  validity  of  a  plane-wave  reflection  coefficient  for  microroughness  than  for 
fine-scale  roughness  in  the  case  of  penetrable  surfaces;  nevertheless,  we  make  this 
assumption  to  avoid  major  complexity  that  will  later  accompany  reflection  in  bistatic 
scattering.  Here  we  apply  it  to  the  monostatic  case  and  isotropic  penetrable  surfaces. 
It  allows  Bragg  scattering  due  to  the  microroughness  that  can  be  written  as 

"V  =  Ro2  r2  cos20  W[2ksin0]  /  n2  ,  (11) 

where  the  power  density  spectrum  is  evaluated  at  the  resonant  wavenumber.  The 
validity  of  Eq.  (11),  however,  is  limited  to  angles  that  are  greater  than  what  might  be 
called  the  "partition  angle",  which  may  be  defined  as  the  angle  subtended  by  a  facet. 
There  are  several  reasons  for  this  limit:  (1)  scattering  would  become  infinite  if  the  angle 
went  to  zero  (i.e.,  the  resonant  wavenumber  would  tend  to  zero  and  the  resonant 
wavelength  and  spectral  density  to  infinity);  and  (2)  the  height  difference  for  widely 
separated  points  on  the  surface  violates  the  small  perturbation  assumption.  Since  the 
facet  size  was  defined  by  g  =1,  it  is  a  reasonable  limit  for  the  validity  of  perturbation 
theory. 


Figure  1  shows  the  contributions  of  mf  and  to  the  scattering  strength.  Note 
that  at  wide  angles  exceeds  mf.  It  also  exceeds  mf  very  near  normal  incidence,  but 
this  part  we  disregard  based  on  the  above  considerations.  We  adopt  the  following 
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scheme  for  adding  mf  and  m  •  m  will  usually  dip  below  mf  for  some  part  of  the 
angular  range;  where  that  difference  (in  dB)  is  the  greatest  (at  0p)  we  start  adding 
to  mf  because  the  abrupt  change  has  the  least  impact.  For  grazing  angles  smaller  than 
this,  the  Bragg  contribution  becomes  increasingly  important,  and  is  added  to  the  rough 
facet  model.  The  result  is  that,  of  the  two  terms,  either  mf  or  m^  dominates  over  most 
of  the  range  with  only  a  small  region  of  transition.  For  very  smooth  surfaces,  where 
there  isn’t  sufficient  microroughness  or  rms  slope  within  the  footprint,  Bragg  may 
dominate  the  rough  facet  model  in  all  but  the  specular  direction.  We  will  see  in  the  next 
section  that  at  wider  angles  the  volume  scattering  term  will  dominate,  and  in  some 
cases  it  will  take  over  before  Bragg  becomes  significant. 

The  Volume  Scattering  ModeJ 

To  estimate  the  contribution  of  mv  coming  from  the  interior  of  the  sediment 
through  the  rough  interface  we  have  adopted  the  Ivakin  and  Lysanov  formalism, 
which  accounts  for  the  volume  contributions  taking  into  account  the  rough  boundary. 
In  their  work  they  assumed  a  uniform  distribution  of  scatterers  within  the  volume  of  the 
sediment.  Since,  to  the  authors  knowledge,  there  is  no  conclusive  evidence  in  the 
literature  that  supports  any  specific  type  of  distribution  for  the  random 
inhomogeneities,  we  have  adopted  the  uniform  distribution  of  isotropic  scatterers, 
which  is  the  simplest  possible  selection  (the  horizontal  stratification  of  the  sediment 
layers  do  not  constitute  the  random  scatters  modeled  in  this  work).  In  the  present 
model,  the  scattering  cross-section  per  unit  volume  (or  volume  scattering  coefficient, 
m  )  is  the  only  free  parameter.  Unfortunately,  most  of  available  datasets  that  can  be 
used  for  testing  the  model  do  not  provide  the  minimum  necessary  geoacoustical 
information  to  apply  the  model  without  requiring  further  "guesstimates".  Hence,  in  all 
the  cases  examined  in  this  work,  other  parameters  beside  m0  need  to  be  estimated. 


In  the  special  case  of  monostatic  backscattering  and  isotropic,  weak,  volume 
scattering  the  Ivakin/Lysanov  model  yields  the  following  expression  for  the  volume 
contribution  to  the  bottom  backscattering  coefficient 


mn  /rV)[n2-sin20’]V2N 

mv  =  2  < - } 

20m  n  cosa 


(12) 


where  p  =  2k 0’=  0-a,  and  a  is  the  local  slope  angle  of  the  interface.  The  angle 
brackets  represent  a  fine-scale  slope  average  over  the  footprint.  T(0)  is  the  local 
amplitude  transmission  coefficient  given  by 

T(0)  =  2  m  cos0  /  [m  cos0  +  (n2  -  sin20)Vz]  . 
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Since  we  want  to  keep  the  algorithm  as  simple  as  possible,  averaging  is  accomplished 
simplv  by  replacing  the  incidence  angle  by  the  incidence  angle  minus  the  rms  slope 
(i.e.  <0’>  =  0-Of).  This  is  a  reasonable  approximation,  since  making  the  effective 

incidence  angle  larger  will  emphasize  scattering  at  high  grazing  angles,  i.e.,  more 
enerqy  will  enter  the  bottom  on  the  upslope  side  than  on  the  downslope  side.  Also, 
this  can  account  heuristically  for  shadowing.  To  further  simplify  the  modeling  process, 
we  use  the  bandlimited  af  value  evaluated  at  6=0  (i.e.,  af0).  Therefore,  the  estimated 
volume  contribution  is  given  by 


m0  ,TVan2-sin2(6-afo)] V* 
2/3m2  n  cosafo 


(13) 


Finally,  we  rearrange  Eq.  (13)  to  produce  the  form 
mV  =  Vy,  *  V(0)  , 


where 

cos40  [1-sin20/n2]Vz  (m+Vn2)4 

V(0)  =  - - -  (15) 

[m  cose  +  (n2  -  sin20)]4 


and 

„  . f  4  (16) 

v  p  cos 6f0  [m  +  Vn  ] 

The  functional  form  V(0)  goes  to  unity  as  0  goes  to  zero  (normal  incidence).  The 
constant  y  is  a  surface  scattering  constant  attributed  to  volume  scattering,  that  is,  it 
is  a  constant  times  an  angular  functional  form  referenced  to  the  surface.  In  addition 
to  showing  the  form  of  the  first  two  terms  in  Eq.  (1),  and  their  combined  effects,  Fig. 
1  also  shows  the  relative  contribution  of  the  volume  term  mv. 

Much  can  be  said  of  the  functional  form  of  V(0),  but  here  we  limit  our  comments 
to  pointing  out  that  it  can  vary  between  cos0  (Lommel/Zeeliger  law  )  at  normal 
incidence  through  cos26  (Lambert’s  law)  at  moderate  grazing  angles,  to  higher  powers 
of  cosine  for  near  grazing.  Many  researchers  have  fitted  backscattering  data  at 
moderate  grazing  with  Lambert’s  law.  In  that  case  becomes  an  empirical  fit  to  the 
theoretically  derived  (Note,  however,  that  /jl  is  not  numerically  equal  because 
V(0)  approaches  normal  incidence  differently  from  Lambert’s  law.)  In  future  work,  pv 
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(and  therefore  the  free  parameter  m0)  will  be  related  to  the  type  of  sediment,  in  a 
manner  similar  to  that  which  is  sometimes  done  for  the  Lambert  constant. 

Since  n  is  complex,  the  calculations  of  R(0)  and  T(0)  are  performed  using  the 
FORTRAN  complex  intrinsic  functions  in  the  computer  implementation.  Allowing  for 
complex  n  constitutes  another  difference  from  the  original  Ivakin/Lysanov  model.  This 
is,  however,  consistent  with  later  work  by  the  same  authors  in  connection  with 
scattering  from  aflat  bottom.19  Complex  n  allows  for  attenuation  in  the  bottom,  which 
becomes  a  key  factor  near  critical  angle  for  a  fast  bottom. 


III.  RESULTS 

The  model  will  now  be  applied  to  a  variety  of  available  datasets.  The  idea  is  to 
examine  the  prospects  for  this  simple  functional  form  to  yield  reasonable  results  when 
fitting  a  wide  range  of  experimental  data.  Table  I  list  the  datasets  chosen  for  this  study, 
the  physiographical  region,  the  surficial  sediment  (as  listed  in  the  original  reference), 
and  the  code  assigned  to  each  dataset.  Most  of  the  data  were  extracted  from  the 
database  developed  by  Li  Zhang.20  To  provide  a  further  distinction  among  the 
datasets,  the  code  name  assigned  to  each  dataset  has  been  modified  from  the 
designation  in  Zhang’s  database. 


Table  I:  Datasets,  locations,  and  codes. 


DATASET 

Buckley  and  Urick24 
Buckley  and  Urick24 
Hines  and  Barry4 
Hines  and  Barry4 
Jin  et  al 21 
Maradino  et  al.26 
Schmidt22 
Schmidt22 

Volovov  &  Zhitkovskiy23 
Zotov  and  Fokin25 


REGION 

Cont.  Rise  (Area  II) 
Cont.  Rise  (Area  III) 
Abyss.  Plain  (Site  A) 
Abyss.  Plain  (Site  C) 
China  Sea 
Abyss.  Plain  (Stn  1) 
Norweign  Sea  (Stn  3) 
Cont.  Slope  (Stn  5) 

N.  Atlant  (Stn  1) 
Clipperton  Fault  (Stn  2) 


SEDIMENT 

CODE 

Silt? 

byuk68  A2 

byuk68_A3 

silt 

hnby92  SA 

silt/sand 

hnby92_SC 

clay 

jnel89 

sand/silt 

mogy85 

scdt69  stn3 

scdt69  stn5 

wzy74 

zvfn85 

Surficial  sediment,  as  listed  in  the  original  reference 
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Table  II:  Geoacoustic  and  experimental  parameters. 


DATASET 

f(kHz) 

lo 

P 

b 

i 

U 

— P- 

a(x10  +  3l 

_b_ 

6f(dea) 

eF(deg) 

byuk68_A2 

2.8 

1545  1585 

1.6 

0.15 

1.2 

1.62 

0.2 

0.04 

0 

byuk68  A3 

2.8 

1530  1510 

1.7 

0.10 

1.2 

1.62 

0.2 

0.04 

0 

hnby92  SA 

0.9 

1545  1680 

1.6 

0.04 

1.2 

1.62 

0.1 

0.13 

0 

hnby92_SA 

1.2 

1545  1680 

1.6 

0.04 

1.2 

1.62 

0.2 

0.04 

0 

hnby92_SA 

2.3 

1545  1680 

1.6 

0.04 

1.2 

1.62 

0.2 

0.05 

0 

hnby92_SC 

2.3 

1545  1550 

1.6 

0.15 

1.2 

1.62 

0.2 

0.05 

0 

jnel89 

2.0 

1530  1510 

1.7 

0.10 

2.0 

1.95 

2.5 

0.06 

5 

mogy85 

0.365 

1545  1650 

2.0 

0.09 

1.2 

1.95 

0.7 

0.33 

0 

scdt69_Stn3 

3.2 

1530  1528 

2.0 

0.10 

2.0 

1.95 

2.7 

0.04 

5 

scdt69_Stn5 

3.2 

1530  1528 

2.0 

0.10 

2.0 

1.95 

2.7 

0.04 

5 

wzy74 

5.0 

1530  1510 

1.7 

0.10 

2.0 

1.62 

0.7 

0.02 

0 

zvfn85 

4.0 

1545  1784 

2.0 

0.70 

2.0 

1.62 

0.9 

0.03 

0 

Table  II  lists  the  parameters  used  for  fitting  the  data  (c0,  c1f  pv  Kp,  a,  and  b), 
and  the  related  surface  parameters  afo,  a  and  0p.  No  precise  geoacoustic  data  are 
available  for  these  datasets.  For  some  of  them  only  a  vague  qualitative  description  can 
be  found  or,  at  most,  some  information  about  the  surficial  sediment  is  provided  only 
for  a  near-by  site.  Therefore,  in  every  case,  we  have  attempted  a  good  fit  by  choosing 
reasonable  values  for  the  input  parameters.  Table  III  lists  the  parameters:  the  Lambert 
constant  pL  (in  dB);  pv  (in  dB),  and  the  resulting  volume  scattering  cross-section  m0; 
as  obtained  by  inverting  Eq.  (8).  In  all  cases,  for  spectral  calculations  of  rms 
roughness,  partition  wavenumber  and  slopes,  the  minimum  wavelength  allowed  in  the 
bottom  roughness  was  set  to  Aa/1000,  where  Xa  is  the  acoustic  wavelength.  The 
largest  wavelength  (footprint)  was  set  to  128m. 

For  a  power-law  spectrum  for  the  ocean  floor  described  by  Eq.  (6),  the 
roughness-related  parameters  needed  in  the  model  are  calculated  as  follows: 

(i)  The  mean  squared  microroughness  (in  m2)  is  given  by  Eq.  (4)  with  =  1 ,  i.e., 


ct^2  =  1  /  4  k2  cos20  .  (17) 

(ii)  The  partition  spatial  frequency  (in  cycles/m)  is  given  by  substituting  o*  from  the 
above  into  Eq.  (7),  and  solving  the  definite  integral,  i.e., 


Kc  =  [KhY  -  y  /  (87ra2COS20)](1/Y)  ,  (18) 
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where  y  =  -2b +  2. 


(iii)  The  mean-squared  slope  of  the  fine-scale  surface  is  given  by  substituting  Kc  from 
the  above  into  Eq.  (8),  i.e., 

6f02  =  [2na2/(r  +  2)]  •  [Kc<*+2>  -  k/*+2>]  .  (19) 


Moderately  Steep  Incidence 


We  first  present 
results  corresponding  to 
data  at  moderate  grazing 
angles  to  normal  incidence. 
Results  for  grazing 
incidence  will  be  examined 
later.  In  all  cases,  along  with 
the  results  of  the  present 
model,  the  prediction 
obtained  by  combining  the 
rough  facet  model  with 
Lambert’s  law  is  shown. 
Both  the  volume  parameter 
and  the  Lambert  constant  pL  are  chosen  to  get  their  respective  best  fit.  In  Fig.  2  the 
fitting  of  dataset  {jnel89}21  is  shown.  This  dataset  corresponds  to  2kHz  in  a  slow 
sediment.  Also  shown  in  this  figure  are  the  relative  contributions  of  the  three 
components  of  Eq.  (1)  to  the  total  scattering  coefficient.  There  seem  to  be  two  well 
defined  regimes:  near  normal  incidence  the  result  is  controlled  by  the  rough  facet 
model  (reflections  at  the  rough  facets)  and  at  moderate  to  large  incidence  angles  it  is 
controlled  by  the  volume  component.  The  only  significant  contribution,  if  present  at  all, 
of  resonant  (Bragg)  scattering  at  the  microroughness  scale  occurs  in  the  transition 
region  (elbow  of  the  curve).  Figures  3a  and  3b  show  results  for  dataset  {scd69}^  at 
3.2kHz  for  two  sites.  In  particular,  in  Fig.  3b,  the  model  mimics  well  the  transition 
between  the  specular  lobe  and  the  middle  angular  range  plateau,  indicating  that  the 
three  components  to  the  scattering  coefficient  are  adequately  combined.  Figure  4 
corresponds  to  dataset  {wzy74}22  at  5  kHz.  The  model  handles  the  transition  region 
quite  well.  It  should  be  noted  that,  for  these  examples  involving  moderate  steep 
angles,  plays  a  crucial  role. 


Table  III:  Fitting  parameters. 


DATASET 

f(kHz) 

n 

Kp_  ±iv(dB) 

ut(dB) 

m, 

byuk68_A2 

2.8 

0.975 

0.15  -33.2 

-25.0 

9.8x1 0'5 

byuk68  A3 

2.8 

1.013 

0.10  -33.7 

-28.0 

6.3x1  O'5 

hnby92_SA 

0.9 

0.919 

0.04  -32.4 

-30.0 

9.3x1  O'6 

hnby92  SA 

1.2 

0.919 

0.04  -37.1 

-35.0 

1.1x1 0'5 

c. 

hnby92  SA 

2.3 

0.919 

0.04  -39.2 

-37.0 

4.9x1 0‘6 

hnby92_SC 

2.3 

0.996 

0.15  -29.0 

-25.0 

2.2x1  O'4 

jnel89 

2.0 

1.013 

0.10  -47.0 

-43.0 

2.2x1 0'6 

mogy85 

0.365 

0.94 

0.09  -33.1 

-30.0 

8.6x1  O'6 

scdt69_Stn3 

3.2 

1.001 

0.10  -28.0 

-27.0 

3.0x1  O'4 

scdt69_Stn5 

3.2 

1.001 

0.10  -34.0 

-33.0 

7.4x1  O'5 

wzy74 

5.0 

1.013 

0.10  -24.1 

-23.0 

1 .0x1 0*3 

zvfn85 

4.0 

0.866 

0.70  -35.2 

-33.0 

4.0x1  O'4 
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Shallow  Grazing  Incidence 


Let  us  now  examine  the  case  of  shallow  grazing  angles.  Figure  5a  shows  results 
for  dataset  {hnby92_SA}4  at  900  Hz.  This  dataset  corresponds  to  a  site  in  the  Sohm 
Abyssal  Plain  which  has  a  fast  sandy  bottom.  The  data  displays  a  hump  near  grazing 
which,  according  to  the  modeling,  occurs  in  the  vicinity  of  the  critical  angle.  Since 
attenuation  is  being  taking  into  account,  total  reflection  does  not  occur.  Lambert’s  law 
is  blind  to  this  anomaly.  Figures  5b  and  5c  show  results  at  1.2  and  2.3kHz  for  the 
same  site.  In  the  three  cases  the  proposed  model  tracks  the  data  fairly  well.  Figure  5d 
shows  results  for  site  C  at  2.3kHz.  In  this  case  there  is  no  hump  and  the  model 
departs  from  data  at  shallow  grazing  angles.  To  "erase"  the  hump  from  the  model,  an 
attenuation  factor  higher  than  the  one  used  for  site  A  was  adopted.  This  choice,  when 
combined  with  the  best  choice  for  ^v,  in  turn  leads  to  a  higher  volume  scattering  cross 
section  than  in  previous  examples.  In  this  particular  dataset,  Lambert’s  law  (also 
shown  in  Fig.  5d)  does  a  better  job  than  the  current  model.  Figures  6a  and  6b  show 
results  for  datasets  {byuk68_A2}24  and  {byuk68_A3}24,  respectively,  corresponding 
to  two  different  sites  with  the  same  frequency  (2.8kHz).  In  the  first  case  the  bottom 
was  modeled  as  a  fast  bottom  to  get  the  hump  at  the  critical  angle.  For  the  second 
case  a  slow  bottom  was  assumed  (see  Table:  II).  In  both  cases  the  fitting  is 
satisfactory.  Lambert’s  law  cannot  track  the  data.  Figure  7  shows  the  fitting  for  dataset 
{zvfn85}2^  at  4kHz  and  Fig.  8  for  dataset  {mogy85}26  at  365Hz.  In  both  cases  the 
model  does  an  excellent  job  in  tracking  the  data,  while  Lambert’s  law  fails  to  mimic  the 
shallow  angles. 

Discussion 


The  role  of  the  different  environmental  parameters  in  fitting  data  varies  with  the 
angular  region.  Near  normal  incidence  the  spectral  parameters  a  and  b  control  the 
lobe  through  6f  and  a  The  impedance  mismatch  plays  a  minor  role  and  the  volume 
contribution  is  negligible.  At  large  incidence  angles  the  fitting  strongly  depends  on  the 
volume  scattering  cross-section  m0  and,  to  a  lesser  degree,  on  both  the  rms  slope  of 
the  large  surface  (it  controls  the  rate  of  decay  near  grazing)  and  the  impedance 
contrast.  For  fast  sediments,  where  a  critical  angle  exists,  the  attenuation  factor  Kp 
controls  the  amplitude  of  the  hump.  In  the  case  of  slow  sediments  the  effect  of  Kp  is 
negligible  for  fitting  purposes.  It  plays  a  role  only  when  backing  out  m0  from  pv 
through  Eq.  (9).  Figures  9  and  10  show  the  sensitivity  of  the  fitting  to  changes  in 
theattenuation  coefficient  Kp  and  in  the  rms  slope  of  the  fine-scale  surface  (6f)  for  large 
incidence  angles  (near  grazing  incidence). 

The  backscattering-derived  volume  scattering  coefficient  (m0)  is  extracted  from 
the  fitting  constant  pv.  Although  the  magnitude  of  the  volume  scattering  coefficient 
(m0)  can  change  by  almost  an  order  of  magnitude  just  by  changing  Kp  for  slow 
sediments,  the  fitting  for  pv  is  almost  insensitive  to  Kp.  On  the  other  hand  for  fast 
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sediments,  the  fitting  is  more  sensitive  to  Kp  and,  therefore,  the  value  of  m0  is  more 
robust.  The  resulting  values  for  m0  are  within  the  order  of  magnitude  of  those  reported 
in  the  literature,  obtained  also,  in  all  cases,  through  indirect  methods.  However,  due 
to  the  crude  averaging  process  adopted  for  the  volume  contribution,  the  value  of  m0 
should  be  considered  a  rough  estimate  of  the  true  value.  Future  work  will  be  oriented 
to  relate  the  free  parameter  (and  therefore  m0)  with  the  type  of  surficial  sediment. 


IV.  CONCLUDING  REMARKS 

A  simple  model  for  estimating  the  bottom  backscattering  strength  for  low  to 
moderately  high  frequencies  has  been  set  forth.  This  model  combines  a 
Helmhoitz/Kirchhoff  two-scale  scattering  model  with  Ivakin  and  Lysanov’s  volume 
scattering  model,  and  with  resonant  scattering  from  the  microroughness.  The  model 
does  not  assume  any  specific  scattering  mechanism  within  the  sediment,  and  aims 
toward  a  rapid  but  reliable  estimation  of  the  scattering  strength  from  a  gross 
classification  of  the  type  of  surficial  sediment  and  physiographic  region  (abyssal  plain, 
continental  rise,  shelf,  etc).  It  contains  only  one  free  parameter  proportional  to  the 
volume  scattering  cross-section  per  unit  volume.  The  model  performs  reasonably  well 
under  a  variety  of  conditions.  The  surface  roughness  seems  to  control  mainly  the 
angular  region  near  normal  incidence  (specular  lobe),  and  the  volume  scattering 
dominates  at  moderate  to  large  angles  of  incidence.  The  resonant  contribution  from 
the  microroughness  dominates  the  transition  angular  range  between  the  two  regions. 
It  will  also  control  the  shallow-grazing  incidence  in  cases  where  the  impedance 
missmatch  prevents  penetration  and  the  volume  contribution  is  weak.  The  attenuation 
in  the  sediment  seems  to  play  a  role  in  scattering  for  fast  sediments.  For  slow 
sediments,  on  the  other  hand,  the  model  is  basically  insensitive  to  variations  in  the 
attenuation. 
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microroughness  starts  to  contribute.  The  scattering 
coefficients,  m^,  mf,  and  mv  are  labelled  by  their 
respective  models,  Bragg,  RFM,  and  volume.  The  overall 
scattering  coefficient,  mbs,  is  given  by  the  solid  curve. 
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